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Abstract. 

As an independent check on previous measurements of the isotopic abundance of oxygen through the Galaxy, we present 
a detailed analysis of the ground state rotational lines of 16 OH and 18 OH in absorption towards the giant molecular cloud 
complex, Sagittarius B2. We have modelled the line shapes to separate the contribution of several galactic clouds along the line 
of sight and calculate 16 OH/ 18 OH ratios for each of these features. The best fitting values are in the range 320-540, consistent 
with the previous measurements in the Galactic Disk but slightly higher than the standard ratio in the Galactic Centre. They do 
not show clear evidence for a gradient in the isotopic ratio with galactocentric distance. The individual 16 OH column densities 
relative to water give ratios of [H 2 O/OH]=0.6-1.2, similar in magnitude to galactic clouds in the sight lines towards W51 and 
W49. A comparison with CH indicates [OH/CH] ratios higher than has been previously observed in diffuse clouds. We estimate 
OH abundances of 10 7 — 10 6 in the line of sight features. 
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1. Introduction IPrantzos et al.lfl99r3) show that the primary/secondary ratios, 

16 0/ n O and 16 0/ 18 0, should fall with decreasing galactocen- 

The OH radical is one of the key oxygen bearing species in the . ,. . , . . , , 

J J , ° r trie distance due to the increased processing rate towards the 

interstellar medium (ISM). Not only is it a good diagnostic of „ , „ ™ , , „ „ ... .. , 

J I — : — a n I -i Galactic Centre. These ratios should also rail with time due to 

physical and chemical conditions (e.g. Goicoechea et al. 2005), u c ., , , . ■ TOl/r 

v J b 1 1 the build up or the secondary elements in the ISM. 

but in addition, it provides a very good way to investigate the 

relative abundance of oxygen isotopes via its isotopologues previouf , measurements of us 0/ i8 in the ISM have mainly 

16 OH, 18 OH and 17 OH. Chemical fractionation reactions that been made yia ^ radio i ines of H?C Q a nd CO towards molec- 

might distort the oxygen isotopic ratios in molecular species uJar douds (see | W ilson & Roodll994T) . However, these lines 

are not thought to be important ganger et al. | | 1984 and this are generally optically thick in the most abundant isotopo . 

means that observations of OH can be used directly to deter- loguef , and go double ratios such as h13 C 16 0/h 12 C 18 are 

mine the ratios 16 0/ 18 and ls O/ 17 0. used TWs rdies Qn accurate knowledge of 12 C/ 13 C which is 

These values are important as they are set by stellar sub j ect t0 chemic al fractionation in molecular species (e.g. 

processing and outflow mechanisms a nd so constrain mod - |L anger et alJl984l) 
els of galactic chemical evolution (e.g. IPrantzos et alJll99r3) . 

ie O is a primary product of stellar nucleosynthesis, pro- A very good way t0 get around these difficulties 

duced directly from the primordial elements H and He (see is t0 use the Q H molecule - several measurements us- 

I Wilson & Matteucci| | l992|) . Both n O and 18 are secondary ing its lg cm A . doubling tr ansitions have been made 

products which require heavier elements from previous nuclear towards the Galact ic Centre llWhiteoak & Gardned flOBlt 

burning for their product ion. Chemical evolution models (e.g. Iwilliams & Gardr^ll98]l although in some circumstances 

* Based on observations with ISO, an ESA project with instruments these can be comphcated by excitation effects within the hy- 

funded by ESA Member States (especially the PI countries: France, perfine levels jBuiarrabal et alJll983h. This problem can be 

Germany, the Netherlands and the United Kingdom) with the partici- avoided by using the far-infrared (FIR) rotational lines which 

pation of IS AS and NASA. should not be affected by hyperfine excitation anomalies. These 
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Fig. 1. Low lying rotational energy levels in the electronic and 
vibrational ground state of the 16 OH molecule, adapted from 
Brow n et al . (1982). The A-doubling of each level is shown but 
the splitting has been exaggerated for clarity. The relative spac- 
ing of the levels in energy is only approximate. The J values for 
each level are shown and the wavelength of each transition is 
given in fim. The transitions used in the analysis presented here 
are shown as thick arrows. 



provide an excellent way to independently check the ratios de- 
termined at radio wavelengths. 

Rotational transitions of OH have previously been studied 
towards the Orion Klei nmann-Low nebula u sing the Kuiper 
Airborne Observatory dMelnicketalJll990t and references 
therein) where both 16 OH and 18 OH were used to con- 
strain the physical conditions o f the source. More recently, 
iGoicoechea & Cernicharol lEoCG) have used the Infrared Space 
Observatory (ISO) satellite to observe OH in the envelope of 
the giant molecular cloud complex Sagittarius B2 (Sgr B2). 
The lowest energy transitions show absorption due to both 
16 OH and 18 OH in Sgr B2 itself and in foreground features in- 
tersecting the line of sight. They indicate that the isotopic ratios 
are broadly similar to the previous radio results. 

In this paper we conduct a more detailed investigation into 
ISO observations towards Sgr B2, with the aim of separat- 
ing the relative abundances of 16 OH and 18 OH in individual 
absorption components in the line of sight. A similar com- 
parison of OH with 17 O H has already been presented by 
IPolehampton et all J2003I) . We have used data from a wide 
spectral survey carri ed out with the ISO Long Wavelength 
Spectrometer (LWS; Clegg et al. 1996) Fabry-Perot mode. A 
combination of prime and non-prime data from the survey al- 
lowed us to increase the signal-to-noise ratio over the standard 
data and derive an accurate and consistent calibration for two 
ground state transitions of 16 OH and one transition of 18 OH. 




-400 -200 200 400 600 800 1000 
^lsr ( km s ~') 



Fig. 2. The 16 OH 2 n l/2 - 2 Tl V 2 .7=3/2-3/2 transition at 53.26 fim 
from non-prime data (black) and prime data (grey). Both 
datasets are binned at 1/4 of their respective resolution ele- 
ments (non-prime: 61 km s _I and prime: 45 km s~'). The gain 
in signal-to-noise in the non-prime data is a factor of ~9. 

After presenting the observations and results, we describe 
a model of the line shape to separate the line of sight absorp- 
tion into 10 velocity components (Sect.^i. In Sect.[5]we assign 
the 16 OH/ 18 OH ratio for each component to a galactocentric 
distance and compare with previous determinations of the iso- 
topic abundances. We discuss the final 16 OH column densities 
in relation to other related species in Sect. |6] 

2. Observations and data reduction 

Sgr B2 was observed as part of a wide spectral survey using 
the ISO LWS Fabry-Perot (FP) mode L03. Unbiased cover- 
age of the whole LWS spectral range (47-196 //m) was car- 
ried out using 36 separate observations with a spectral res- 
olution of 3CM-0 km s -1 . No other object outside of the 
Solar System was observed over the complete LWS spectral 
range in this way. Results from this survey have been pre - 
sented bvlCeccarelli et alJ ll2002h: IPolehampton et alJ l|2Q02a); 
IVastel et alJ d2002l) : IPolehampton et all (120031 120051) . Sgr B2 
was also extensively observed by the LWS FP in narrow wave - 
length scans using the L04 mode (see lGoicoechea et alJ2004l) . 
We have used several of these L04 observations to improve the 
signal-to-noise ratio in the L03 data for I8 OH. The ISO TDT 
numbers for all the observations used in this paper are detailed 
in TableEO 

The LWS beam had an effective diameter of approximately 
80" JOrvet alJ2 003) and L03 observations were pointed at co- 
ordinates a = 17 h 47 m 21.75 s , 6 = -28°23'14.1" (J2000). This 
gave the beam centre an offset of 21.5" from the main FIR 
peak at Sgr B2 (M) - this pointing was used to exclude the 
source Sgr B2 (N) from the beam. The additional L04 obser- 
vations were pointed directly towards the nominal position of 
Sgr B2 (M), but careful comparison showed no significant dif- 
ference to the L03 observations. Each observation had a spec- 
tral sampling interval of 1/4 resolution element with each point 
repeated 3-4- times in L03 mode and 1 1-15 times in L04 mode. 
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Fig. 3. Data for the I6 OH cross ladder ( 2 Tli /2- 2 ^3 /i) transitions, 7=1/2—3/2 at 79 fim (left; binned at 1/4 resolution element) and 
J -3/2-3/2 at 53 fim (right; binned at ~ 1/6 resolution element), plotted as line to continuum ratio. The errors shown represent the 
combined statistical and systematic uncertainty in each bin. The two A-doublet components have been averaged for the 53 fim 
line. The best fitting model accounting for the absorbing clouds in the line of sight is shown as the solid line (see text). 



In addition to the primary spectral survey data, each line 
was generally included in at least one other L03 observation. 
This occurred because the LWS used 10 detectors, which al- 
ways recorded data in their own spectral ranges (the prime 
observations consist of data from a single detector for which 
instrument settings were optimised). We have included these 
'non-prime' data to further improve the signal-to-noise ratio. 

Each observation was carefully reduced using the LWS 
offline pipeline (OLP) version 8 (for FP data the difference 
between OLP version 8 and the latest version 10 is not sig- 
nificant). Further processing was then carried out interac- 
tively using routines that appeared in the LWS In teractive 
Analysis package version 10 (LIA10: lLimet al.l2002l) and the 
ISO Spectral Analysis Package (ISAP: ISturm et alll998h . The 
method included determination of accurate dark currents (in- 
cluding stray light), interactive division of each mini-scan by 
the LWS grating response profile s and careful scan by scan re- 
moval of glitches (seelPolehampt orTet alJ2002all2003l) . 

Before co-adding the data for each line, the wavelength 
scale of each observation was corrected to the local standard 
of rest (LSR). Non-prime observations were carefully checked 
against the equivalent prime measurement to align the line cen- 
tres. This is particularly important for non-prime data measured 
with the LWS long wavelength FP (FPL) below 70 fim (out- 
side its nominal operating range) because the standard FPL 
wavelength ca libration did not include short wavelength data 
JOrvetaljf2 003). For the 16 OH J -3/2-3/2 cross ladder transi- 
tion at 53 fim (Fig.^shows the energy levels for OH), we have 
used only non-prime FPL data. This is because it has a sig- 
nificantly higher signal-to-noise ratio than the prime data (due 
to the higher tr ansmission of FPL compar ed to the short wave- 
length FP, FPS; Polehampton et alJ2002bl) . The absolute wave- 
length alignment was difficult to determine purely by compari- 
son with the noisy prime data and so a velocity shift was al- 
lowed as a free parameter in the modelling (see Sect. 14. 1> . 
Figure|2]shows the resulting shift cross-checked with the prime 



data. The gain in signal-to-noise achieved by using FPL rather 
than FPS is approximately a factor of 9. The only cost in using 
these data is a reduction in spectral resolution from 45 km s _1 
(for FPS) to 61 km s 1 (for FPL). 

After co-addition, the continuum around each line was fit- 
ted with a 3rd order polynomial baseline which was then di- 
vided into the data to obtain the relative depth of the lines below 
the continuum. This effectively bypassed the large systematic 
uncertainties in the multiplicative c alibration factors need ed to 
obtain the absolute flux scale (see ISwinvard et aDll9 98). The 
remaining errors are due to detector noise, uncertainty in the 
dark current and the polynomial fit. 

At the resolution of the LWS, two components are visible 
for each OH transition due to the A-doublet type splitting of 
each rotational level. However, further hyperfine splitting is not 
resolved. The two A-doublet components showed good agree- 
ment in the data, and so where they were well separated (the 
7=3/2-3/2 cross ladder transition at 53 fim for 16 OH and the 
7=5/2-3/2 transition at 120 fim for 18 OH), they were co-added 
to further increase the signal-to-noise ratio. 

3. Results 

Figure ^ shows the low-lying rotational transitions of 16 OH, 
most of which are included in the spectral survey range. 
The strongest lines observed in the survey were due to 
the 2 Ylj,/2 7=5/2-3/2 transition from the ground state at 
119 fim. These lines show almost complete absorption of 
the FIR continuum in the range -150 to +100 km s _1 (see 
Goic oechea & Cernicharoll2002t) . Due to the large depth of 
the lines, the shape was strongly affected by the transient re- 
spon se of the LW S detectors. These detector memory effects 
(see lLlovdbOol) meant that successive repeated scans underes- 
timated the depth and the continuum level following each line. 
Therefore, the lines were not used in the analysis presented 
here. There are two remaining transitions from the ground rota- 
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Fig. 4. Co-added data from the two A-doublet components of 
the 18 OH 2 Il3/2 7 — 5/2-3/2 transition. The data were binned 
at 1/4 resolution element and the errors represent the combined 
statistical and systematic uncertainty in each bin. The best fit- 
ting model accounting for absorbing clouds in the line of sight 
is shown (see text). 



tional state (see Fig.Q) that occur between the 2 rLj/2 and 2 I\-\j2 
ladders: 7=1/2-3/2 at 79 /mi and 7=3/2-3/2 at 53 /mi. These 
lines are also broad but with much lower optical depth. They 
are shown in Fig. [3] The broad absorption is due to features 
between the Sun and Galactic Centre associate d with galactic 
spiral arms that cross the line of sight (e.g. lGreaves & Williams! 
11994 . 

Several I6 OH transitions between higher energy rotational 
levels are observed in the survey at the velocity of Sgr B2 it- 
self but do not occur in the line of sight clouds. These lines 
origi nate in the envelope o f Sgr B2 and have been modelled by 
IGoicoechea & Cernichard J2002I) . They show that the excited 
OH originates in clumpy photodissociation regions (PDRs) on 
the edge of the Sgr B2 complex at temperatures of 40-600 K. 
In this paper, we concentrate on the broad absorption observed 
in the ground state lines at 79 and 53 /mi. 

Due to the lower abundance of I8 OH, the fundamental 
ground state 2 fl3/2 7=5/2-3/2 transition is much weaker and the 
cross-ladder and higher energy lines are not detected. Figure^] 
shows the 7=5/2-3/2 line after co-adding the two resolved A- 
doublet componen ts. These observa tions were previously pre- 
sented by Polehampt on et al . (2003). Here we have used addi- 
tional data from the L04 mode, which showed good agreement 
with the L03 data. 

4. Modelling of OH in the line of sight 

4.1. High resolution model 

At the spectral resolution of the LWS, the line of sight compo- 
nents are blended together into a single broad absorption. We 
have modelled the line shape using higher spectral resolution 
measurements to fix the velocities and widths of each compo- 
nent. Comparison of high resolution spectra tracing molecu- 
lar and atomic species in these clouds show very similar line 



widths and velocities (C O; e.g. lVastel etaD(l2002i). H?CO; e. g. 
Mehringer et all lll995l) . HI; e.g. Idarwood & Dickevl Jl989k 
These tracers also show a very similar velocity structure to 
observati ons of the A -doubling lines of OH at cm wave- 
lengths llMcGee et alJll97(ilWhiteoak & Gardner! 1 9751 1 1 98 it 
IWilliams & Gardnerll98lUBuiarrabal et alJl983l) . 

We have assumed that the H i measurements are a represen- 
tative tracer_ajidba^ed_our model on the published parameters 
from iGarwood & Dickevl Jl989l) (after correcting an error in 
their Table 2: 11 km s _1 should read 1.1 km s _1 ). This method 
has alrea dy been used to suc cessfully fit the absorption of CH 
and CH2 llPolehampton et all 2005) an d is s i milar to the method 
used to fit 1 and C 11 lines by Va stel et al.l (|2002). 

The Hi observations were made using the VLA pointed 
close to Sgr B2 (M). In our model, the Hi data were used to 
fix the velocity and line width of 10 line of sight features with 
each component assumed to have a Gaussian line shape. The 
optical depth at each velocity, t, was then adjusted as a free 
parameter and the line-to-continuum ratio calculated from, 



/ = l c exp (-t) 



(1) 



where l c is the intensity of the continuum. The resulting spec- 
trum was then convolved to the resolution of the LWS. In order 
to put the best constraint on the line shape, the two 16 OH lines 
were fitted simultaneously with the 18 OH line. Optical depths 
were varied in the model for the 79 /mi line and then used to 
calculate the shape of 53 fim line assuming that both lines trace 
the same column density. This must be true as they originate in 
the same lower level. Column densities were calculated assum- 
ing a Doppler line profile with Maxwellian velocity distribution 
re.g.lSnitzeJl97^. 



to Av 



2VhT2 AfjASjgi/gj 



(2) 



where Nj is the column density in the lower level in cnT 2 , To is 
the optical depth at line centre, Av is the line width in km s _1 , 
A t j is the Einstein coefficient for spontaneous emission in s _1 , 
Aij is the wavelength in /mi and g t is the statistical weight of 
state i. The line wavelengths and Einstein coefficients used (av- 
eraged over the unresolved hyperfine structure) are shown in 
TablelA2l 

Further free parameters were added to the fit describing 
the 16 OH/ I8 OH column density ratio, allowing the shape of 
the 120 /mi I8 OH line to be calculated. Clouds which are 
thought to reside at similar galactocentric distances were forced 
to have identical isotope ratios - this affects the features at 
-108/ - 82 km s _1 , -52/ - 44 km s" 1 and +53/ + 67 km s -1 . 
In order to account for any drift in the wavelength calibration 
between the lines, a free velocity shift was allowed for each 
line, giving a total of 20 free parameters in the fit. The most 
important factor constraining the velocity shift was the rela- 
tively deep and narrow absorption at the velocity of Sgr B2 in 
the 18 OH line. The Numerical Recipes multi-parameter fitting 
routine, 'amoeba' (Pres s et al J 19921) . was used to minimise^- 2 . 

In order to ensure that the algorithm did not converge at 
a false minimum, the initial conditions were set close to their 
best values and the routine was restarted at the first convergence 
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Table 1. Column densities and 16 OH/ I8 OH ratios for the best fit car ried out to the 16 OH and 18 O H ground state lines. Columns 1 
and 2 show the Gaussian fit parameters from the H i observations of lGarwood & DickevMl989l) . 



Velocity 




FWHM 




jv( 16 oh) 


16 OH/ 18 OH 


(km s" 1 ) 




(kms" 1 ) 




(10 15 cm" 2 ) 




-108 ±2 




7±4 




0.56 ± 0.39 


460 +1J0 


-82 ±4 




28 ± 12 




9.0 ±3.9 


460+ 90 

too 50 


-44.0 ± 0.6/(-52 


±6) D 


8 ± 2 (17 ± 


5) 


3.5 ± 1.8 


450+ ioo 


-24.4 ± 0.8 




14 ±2 




5.4 ± 1.8 


370+^° 


+ 1.1 ±0.4 




19 ± 1 




7.6 ±3.4 


4 3 0+ 80 


+ 15.7 ±0.3 




7± 1 




3.6 ±2.3 


360 +8 ° 


+31 ±4 




21 ± 14 




5.4 ±3.8 


540+ 100 


+66.7±0.5/(+52.8 


± 0.9)* 


16± 1 (11 ± 


:2) 


32!* 


320+70 

-30 


Total 








67.0 





" These two components are too close to separate in our model and the best fit line shape requires only the -44.0 km s 1 component. 

* The Hi data resolve 2 components in Sgr B2. However, our fit requires only one of these (at 66.7 km s~') to reproduce the ISO spectrum. 



point. Solutions with negative absorption were avoided by set- 
ting^ 2 to be high when the optical depth went below zero. 

The best fitting model is shown plotted with the data for 
each line in Figs. [3]and|4] The final 16 OH column densities 
and 16 OH/ 18 OH ratios are shown in Tabled In the line of sight 
clou ds where no absorption is observ ed from higher energy lev- 
els jGoicoechea & Cernicharol2002l) . the ground state popula- 
tion is a good measure of the total OH column density. The 
fitted components are all optically thin except at the velocity of 
Sgr B2 itself, where optical depths of 2.5-3.3 were found for 
the 16 OH lines. However, these are not high enough to break the 
assumption of purely Doppler line profiles used in the model 
and calculation of column densities. 

The best fitting velocity shift applied to the 16 OH 7=1/2- 
3/2 line at 79 fim was 8.8 km s" 1 and for the 18 OH 7=5/2- 
3/2 line at 120 jum was -1.9 km s . These shifts centred 
the deepest absorption on the Hi component at +67 km s _1 
and are withi n the unc ertainty in FP wavelength calibration 
(< 1 1 km s-': lGrvet all2003l) . However, for the 16 OH 7=3/2- 
3/2 line at 53 fim using data taken with non-prime detectors, a 
larger shift of 28 km s _1 was required. This shift is due to the 
fact that this line was observed with FPL outside of its nominal 
range (see Sect. |2j. In order to determine if this large velocity 
offset was reasonable, we compared the shifted data to the nois- 
ier prime observation (in which the accurac y of the wavel ength 
calibration should be better than 6 km s ^ iGrvet alJ200l .The 
two line components show good agreement (see Fig.[5J and the 
line shift is consistent with that f ound for other spec tral lines 
observed using FPL below 70 fim (Polehampton 2002). 

The final fit shows that only 8 of the 10 velocity compo- 
nents are necessary to reproduce the observed line shape. At 
the velocity associated with Sgr B2, the width of the absorp- 
tion in the 18 OH line is too narrow to allow strong absorption 
by both the +53 and +67 km s components observed for H i. 
This is consistent with t he results obtained by fitting the ground 
state line of CH (Polehampt on et alJ2005l) . Also, the two com- 
ponents at -52 and -44 km s _1 are too closely spaced to be 
separated in our fit and significant optical depth was found only 



at -44 km s 1 (in Tabled only one value centred at -44 km s 1 
is given for the two components). 

We have also investigated the effect of further reducing the 
number of fit components in order to increase confidence in 
the final column densities and ratios. The minimum number of 
components th at can reasonably b e fitted to the spectrum is 4 
(e.g. following |Niieidetalj[2000 ) . We re-ran our model us- 
ing 4 velocity components corresponding to the ranges used by 
Neufeld et al. and widths estimated from their H^O spectrum 
observed with the Submillimeter Wave Astronomy Satellite 
(SWAS). These 4 empirically determined components can 
broadly reproduce the observed line shapes. The fit gives 16 OH 
column densities and 16 OH/ 18 OH ratios that are within the er- 
rors of the results from Tabled summed in the relevant veloc- 
ity ranges. The 16 OH/ 18 OH ratios obtained are, 480 (centred 
at -80 km s -1 ), 470 (centred at -40 km s -1 ), 510 (centred at 
+ 10 km s~') and 290 (centred at +65 km s _1 ). However, we 
have used the high resolution H i observations as a basis for our 
model because this can not only fully describe the line shapes 
but also allows the contribution of clouds at different galacto- 
centric distances along the line of sight to be disentangled. 

4.2. Errors on fitted parameters 

The number of closely spaced velocity components in the fit 
made the modelling difficult as their separation was less than 
the resolution of the LWS. This meant that variation in one 
component could be compensated for by changing another, re- 
sulting in a relatively large uncertainty in each fitted optical 
depth. In addition, the results for each component are likely to 
be an average over several closely spaced narrow features such 
as those obs erved showing CS ab sorption with velocity widths 
~1 km s _1 JGreaves & Williamsl fr994). However, this would 
only alter the total column densities if a few of the narrow fea- 
tures had very much higher optical depths than the others, and 
this does not appear to be the case in the CS data. 

In order to estimate the uncertainty on each fitted param- 
eter we performed a Monte-Carlo error analysis. The best fit- 
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ting model determined by minimising x 2 was used to generate 
a set of synthetic spectra where each point had a mean value 
equal to the best fit and standard deviation equal to the origi- 
nal data error. We re-fitted each synthetic spectrum using the 
original fitting method and analysed the resulting dataset for 
each parameter. The results show that the errors on neighbour- 
ing optical depths are strongly correlated, but that these are not 
correlated with the uncertainty in 16 OH/ 18 OH ratios (except for 
the Sgr B2 component). This is due to the fact that the best 
fit isotopologue ratios depend more on the overall line shape 
than on the relative optical depth in neighbouring components. 
The final uncertainty in column density is shown in Table ^ 
as a combination of lcr e rrors in H i line width (as quoted by 
iGarwood & Dickevlll989l) and the modelling errors in optical 
depth. The errors in 16 OH/ 18 OH ratio were determined directly 
from the Monte-Carlo analysis. 



4.3. Comparison with previous results 

The column densit y at the velocity of Sgr B2 has previously 
been calculated bv lGoicoechea & Cernicharol lEoCG) using ISO 
FP data observed in the L04 mode. They used a radiative trans- 
fer model and accounted for the populations in higher energy 
levels up to 420 K above ground to determine A^( 16 OH) = 
(1.5-2.5) x 10 16 cirT 2 . Lower resolution observations with 
the LWS grating mode show that OH is widespread across 
the whole Sgr B2 regi on with column de nsities in the range 
(2-5)xl0 16 cur 2 llGoicoechea et al.l l2004). Our value associ- 
ated with Sgr B2 from Table □ gives N( 16 OH) = (3.2+°£) x 
10 16 crrT 2 in the ground state level. This is slightly higher 
than the previous FP result but fits into the general picture 
reasonably well. We also find a slightly higher value than 
iGoicoechea & Cernicharol ^002) for JV( 18 OH): (9.8 ± 1.6) x 
10 13 cirT 2 compared to (6 ± 2) x 10 13 cirT 2 . 

The column density of 16 OH has also been measure d from 
observations of its A-doubling lines at 18 cm. iBieg in g l J 1974 
used the interferometer of the Owens Valley Radio Observatory 
(OVRO) to calculate an integrated column density over the ve- 
locity range 40-89 km s _1 equal to 1.2 x 10 17 cirT 2 . This would 
indicate that the FIR lines underestimate the column density 
by a factor of 4. However, if the excitation temperature of the 
18 cm lines is less than the assumed 20 K (see the discussion in 
IStacev et al.ll987l) . the radio column density could be an over- 
estimate. Also, the 18 cm lines may sample a different com- 
ponent of the envelope as the depth at which the continuum is 
emitted occurs deeper into the cloud than in the FIR region, 
and t he OVRO synt hesised beam (3.25') is larger than that of 
ISO. iBieginghl 97(h estimates that the column density in the 
negative velocity features is smaller than the positive velocity 
values by about a factor of four. 



The total column density of 18 OH over the whole line of 
sight in our fit was (1.8 ± 0.2) x 10 14 cirT 2 and this com- 
pares favourably with p revious observation s using the Kuiper 
Airborne Observatory: iLugten et alJ Jl986h found a total col- 
umn density of A?( 18 OH) > 2 x 10 14 cm' 2 . 



o 



o 



800 

600 
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2 4 6 8 10 

Distance from Galactic Centre (kpc) 

Fig. 5. Values of 16 OH/ 18 OH from Table [T]plotted against dis- 
tance from Galactic Centre (large circles with cross). The ve- 
locity of each component is indicated in km s" 1 next to each 
point. The compilation of values from iKahand i 1 994h is plot- 
ted for c omparison, as well as th e best fitting galactic gradient 
found by Wilson & Rood! d 1994 . 



5. 16 OH/ 18 OH ratio 

Table^gives the best fitting I6 OH/ 18 OH ratios. In order to com- 
pare with previous measurements of the isotopic ratios, a galac- 
tocentric distance must be assigned to each feature. Several of 
the velocity components have been well established to come 
from galactic spiral arm and Galactic Centre features such a s 
the Galactic Bar at large negative velocities JScoville| | l972|) . 
the 3 kpc expanding arm at -44 km s" 1 (e.g. lBurkelll9 65). 
the 4-5 kpc arm at -24 km s _1 JMenon & Ciottilll970l) and 
Sgr B 2 (M) itself at +50-70 km s~' (e.g. lMartm-Pintado et alJ 
1990p. A ssociations for the remainin g features have been pro- 
posed by|Greaves_& Williams ( 1994). In this case, the 1, 16 and 
31 km s _1 features would be associated with local spiral arms. 
However, it has been suggested that there may be a contribution 
at k m s _1 by gas associated with Sgr B2 (e.g lGardner et alJ 
1988). 

Figure |5] shows our results plotted ag ainst their distances 
from the Galactic Centre taken from iGreaves & Willi ams 
(1994). Compilations of many observations of H2CO through 
the Galaxy indicate that there is a gradient in the 16 Q/ 18 Q ra- 
tio with decreasing values towards the Ga lactic Centre ifrbsil 
ll982tlWilson &Roodll994tlKahanell994 . The best fit gradi- 
ent found bv lwlkOTi"&RoodT(ri994l) is overplotted in Fig. |5]as 
well a s the compilation of H2CO and CO results from Kahane 
dl994l) . 

The three most uncertain distances are for the local com- 
ponents. At 31 km s we find a relatively high value of 
I6 OH/ 18 OH, consistent with its location close to the Sun. At 
1 km s _1 , our value is in goo d agreement with that found from 
OH A-doublet lines of 440 (Bui arrabal et"aT1l 1983b . However, 
both the 1 and 16 km s _1 components have rather lower val- 
ues than might be ex pected for local gas. The gradient of 
IWilson & Rood! {l994) predicts a local value of 560 ± 25, al- 
though it has been measured to be even higher using HCO + 
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absorption by Lucas & Lisztl dl998l) . giving 672 + 110. This 
may indicate that the local absorbing gas towards Sgr B2 has 
atypical isotopic abundances. An anomalousl y low 1 2 C/ 13 C ra- 
tio has also been found in this gas bv lGreavesI d 1 995 ) . She rules 
out an error in the distance estimate and attributes the extra 
13 C to enrichment by stellar ejecta containing secondary nucle- 
osynthesis products. In this case it would be natural to also find 
a higher abundance of 18 OH. 

In comparison, the standard Solar System abundance (de- 
rived from Standard Mean Ocean Water) is ( i6 0/ 18 0)smow = 
498.7 + 0.1 teaertschilll973) . The value in the outer layers 
of the Sun has be en measured to be 16 0/ 18 Q = 440 + 50 
llHarris et all 19871) and in the Solar wind, Collier et all dl998h 
found 16 0/ 18 = 450 + 130. This should represent the con- 
ditions when the Sun was formed and models predict that it 
should be higher than the current ratio in the local ISM due 
to th e build up of seconda ry elements from stellar processing 
(e.g. IPrantzos et aljfl996l) . In order to explain the low Solar 
System abundances compared to the ISM, it has been pro- 
posed that the Sun was formed closer to the Galactic Centre 
dWielen & Wilson! 19971) . 

In t he Galactic C entre, the isotope ratio is generally taken to 
be 250 dWilson & Roodll994l) . This low value has been derived 
from measurements of the A-dqublet lines of 16 OH and 18 OH 
(e.g. IWiUiams & Gardnerl|l98ll IWhiteoak & Gardnerlfl98ll) . 
However. iBuiarrabalet all dl983l) show that the 16 OH/ 18 OH 



opacity ratio derived from the radio lines is only a good mea- 
sure of the actual abundance ratio in cold clouds with kinetic 
temperatures <20 K. This is due to excitation anomalies in the 
hyperfine levels caused by rotational pumping by the FIR lines. 
This certainly occurs at the velocity of Sgr B2 as can be seen 
from the higher energy lines of 16 OH which show absorption 
with in the 2 Hmj ladder but emis sion within the 2 Yli/2 ladder 
( Goi coechea & Cer nicharo 2002). This would cause the previ- 
ously derived isotopic ratios to underestimate the true value. At 
the velocity of Sgr B2, our fit giv es a ratio of 320 + 3i. Thi s ap- 
pears to confirm the suggestion of Buiarra bal et alJ dl983l) that 
the ratio in the Galactic Centre has been underestimated (as- 
suming the same excitation conditions exist for 16 OH and I8 OH 
and the ground state ratio accurately reflects the isotopic abun- 
dance). We also find a high value for the component attributed 
to the Galactic Bar (-108 km s _1 ) of 460 + ?°. Buiarrabal et ail 
d 1983b find a ratio of ^370 in this component. 

In contrast, measurements of other molecules in Sgr B2 
seem to agree with the lower radio OH values. However, none 
of these methods are easy to directly interpret and so may also 
underestimate the true value. Both H2CO and CO must be used 
in a double r atio, r elying on an accurate value of 12 C/ 13 C. 
iHenkel et all dl983h calculate 16 O/ 18 O~200 from measure- 
ments of H 2 3 C 16 0/H 2 2 C 18 0, although the final value is uncer- 
tain due to the correction for radiation trapping in H 2 2 C lb O 
(used to determine the associated 12 C/ 13 C ratio). Observations 
of CO ind icate a value of 250±100 in the Galactic Centre 
dPenziasll98ll) but this m a y be affected by chem ical fractiona- 
tion of 13 C dPenziasIl 9*831) . Idardner et aljfl989T) have observed 
CH^ 6 OH and CH' 8 OH towards Sgr B2 leading to an isotopic 
ratio of 210+40. This does not include a cor rection for radia- 
tion trapping which would increase the ratio. iNummelin et alJ 



Table 2. Ratio of our derived 16 OH col umn densities with those 
of CH from Polehampt on et all d2005l) for each of the velocity 
components from Tabled 



Velocity 


[ 16 OH/CH] 


(km s" 1 ) 




-108 


6.2 ± 6.5 


-82 


35 ±21 


-52/-44.0 


19 ±15 


-24.4 


28 ±11 


+ 1.1 


38 ± 18 


+ 15.7 


12 ±8 


+31 


36 ±35 


+52.8/+66.7 


34 ±7 



d2000l) calculate even lower ratios from S 16 0/S 18 (120^ 8 ) 



and S lb 2 /S lb 18 (112+ 2b ). They explain these by assum- 
ing that the observed emission originates in the central regions 
of the Sgr B2 complex where ls O could be enhanced by the 
ejecta from massive stars. Our FIR absorption measurements 
trace gas in the outer parts of the Sgr B2 envelope. 

Without low values at the Galactic Centre, our results show 
a much less convincing gradient than that obtained from the 
CO and H2CO measurements. However, we have low num- 
ber statistics and the distances of some components from the 
Galactic Centre are not very well known. There also must be 
some effect from a real variation in clouds at the same galac- 
tocentric distance, as indicated by the anomalous 16 OH/ 18 OH 
and 12 C/ 13 C at km s . This means that a larger sample is 
necessary to accurately confirm the trend through the Galaxy. 

6. Comparison with other species 

Several other related species have been observed in the line of 
sight clouds towards Sgr B2. Observations of absorption due 
to the H 2 b O and H ^ 8 Q ground state ro tational lines have been 
made using SWAS (INeufeld et all2000l) . Although H' b O com- 
pletely absorbs the continuum, they calculate column densities 
in three line of sight velocity ranges from H 2 s O. Both OH and 
H2O are formed by the dissociative recombination of H30 + 
(H30 + is also observed to show absorpt i on in the line of sight 
clouds: iGoicoechea & Cernicharol200ll) . INeufeld et alJ §002) 
show that if recombination of HsO + is the only production 
mechanism, the ratio of OH to H2O gives good constraints on 
the branching ratios. However, they also show that in warm gas, 
neutral-neutral reactions can alter the balance. A further effect 
is the inclusion of grain surface reactions, which can increase 
the gas -phase abundance of H2O whilst OH is relatively unaf- 
fected dPlume et all20 04). 

Us ing the orfho-H 2 8 column densities from lNeufeld et all 
(2000) with an ortho-to-para ratio of 3 and our le O/ ls O ratios 
averaged in their velocity ranges gives [H2O/OH] ratios of 0.6- 
1.2 towards Sgr B2. These values are a similar order of magni- 
tude to those previously d erived for line of sigh t absorption fea- 
tures t owards W51 (-0.3 ; INeufeld et al.l2002l) and W49 (0.26, 
0.43; IPlume et al.l 2004). indicating a high branching ratio for 
H2O and/or the necessity of including other processes such as 
gas-grain interactions. In gas at the velocity of Sgr B2 itself, 
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iGoicoechea & Cernichard J2002) have derived even higher val- 
ues of [H2O/OH]~l-10. In this case, there is probably a warm 
gas component present and neutral-neutral reactions could be 
important. 

In diffuse clouds, CH has been observed to be tightly corre- 
lated with OH T with a column density ratio, [OH/CH]~3 (e.g. 
iLiszt & LucaslEoO^ . Absorption from the ground state rota- 
tional transition of CH was observed as part of the ISO spec- 
tral survey towards Sgr B2 and has been fit ted in a similar 
way to the OH lines (Polehampt on et alJl2005l) . Table|2]shows 
the ratio of our fitted I6 OH column densities those of CH in 
each velocity component. These have a large uncertainty but 
seem to indicate values that could be up to ten times higher 
than that seen in diffuse clouds. This increased ratio cannot be 
caused by CH because the abundanc es found are v ery simi- 
lar to diffuse cloud models (Poleham pton et alJ l2005). We can 
estimate the abundance of OH in the line of sight clouds us- 
ing molecular hydroge n column densities in the range 4 to 
14xl0 21 cirT 2 (see ldreaves & Nvmanll996l) . This leads to OH 
abundanc es 10~~ 7 -10~ 6 in these feature s. At the velocity of 
Sgr B2. lGoicoechea & Cernicharol ( l2002l) find an abundance of 
(2-5)xl0~ 6 . These abundances are much higher than found in 
diff use clouds (mod els of diffuse clouds predict X(OH)~ 1CT 8 ; 
Ivan Dishoeck & Blackl 19861) and lead to the high [OH/CH] ra- 
tio. Such high abundances of OH can be found in PDR regions. 

7. Summary 

In this paper we have presented an analysis of 16 OH and 18 OH 
FIR rotational lines towards Sgr B2. Absorption at different 
velocities along the line of sight allowed us to calculate the 
16 OH/ I8 OH ratio between the Sun and Galactic Centre. 

The advantages of this technique for determining the iso- 
topic ratio of oxygen are: 

- Chemical fractionation effects are not important for 16 OH 
and I8 OH. Therefore the molecular column density ratio 
directly measures the isotope abundance. 

- There are optically thin lines present from both species. 

- The final ratio does not depend on the value of I2 C/ I3 C as 
for previous H2CO and CO measurements. 

- Uncertainties in the population of hyperfine states do not 
affect the final ratios (as they do for radio observations of 
A-doublet transitions). 

- Only ground state lines are seen for the line of sight clouds, 
showing that excitation effects are not important in these 
clouds. 

- All observations come from a single dataset with consistent 
and stable calibration over lines of both isotopologues. 

The main disadvantage is the relatively low spectral resolu- 
tion compared to the velocity separation of line of sight compo- 
nents. This introduced errors for individual components in the 
line shape modelling. However, the data used here allowed us 
to maximise confidence in the observed line shapes by combin- 
ing several observations and with the good coverage of base- 
lines around the lines. 

We find I6 OH/ 18 OH ratios broadly consistent with previous 
isotopic abundances, although our results do not provide such 



clear evidence for a gradient of the ratio through the Galaxy 
(but they do not rule it out). In velocity components associ- 
ated with the Galactic Centre, we find slightly higher ratios 
than previous results. This could be due to an underestimate 
in the previous radio observations of OH caused by pump- 
ing of hyperfine level s by the FIR rotational transitions (see 
Buiarraba l et al.ll983h . 

In comparing the total 16 OH column densities with those 
of H2O towards Sgr B2, we find similar magnitude abundance 
ratios as have been seen in other lines o f sight towards W5 1 and 
W49 jNeufeld et al. 2002; Plu me et alJ20 04). A comparison of 
I6 OH with the column density of CH shows a wide scatter and 
large error, but appears to indicate higher ratios that previously 
observed in diffuse clouds. 

This study shows that FIR rotational lines are an extremely 
useful tool for examining the variation of isotopic ratios in the 
ISM. The column densities of OH can be accurately derived 
as the method avoids the anomalous excitation that may be 
present in radio A-doublet lines. Ideally, future observations 
at higher spectral resolution are required to improve on the 
results. This may be possible using future telescopes such as 
SOFIA. 
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Appendix A: Observations and wavelengths 
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Table A.l. Log of the observations used. 



Transition 


Observing 


ISO 


LWS 


Spectral 


Repeated 




Mode 


TDT 
Number 


Detector 


Resolution 

(kms- 1 ) 


scans 


16 OH 












2 n 3/2 7=5/2-3/2 (119 fim) 


L03 


50601013 


LW2 


33 


3 


2 n 1 / 2 - 2 n 3 /2 


L03 


50400823 


SW2 


61 


3 


7=3/2-3/2 (53 fim) 


L03 


50601013 


SW2 


61 


3 




L03 


50700707 


SW2 


61 


3 




L03 


50800317 


SW2 


61 


3 


2 n 1 /2- 2 n 3 /2 


L03 


50900521 


SW5 


40 


3 


7 =1/2-3/2 (79 yum) 


L03 


50600814 


SW5 


40 


3 




L03 


50600603 


SW4 


40 


3 




L03 


50700208 


SW4 


40 


3 




L03 


83800606 


SW4 


40 


4 




L03 


50600506 


SW4 


40 


3 


ls OH 












2 n 3/2 7=5/2-3/2 (120 ixm) 


L03 


50601013 


LW2 


33 


3 




L03 


50700610 


LW2 


33 


3 




L04 


49401705 


LW2 


33 


15 




L04 


46201118 


LW2 


33 


11 




L04 


46900332 


LW2 


33 


14 



